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G\ ; Abstract 

We derive simple and physically transparent expressions for the con- 
tribution of the strong interaction to one nucleon removal processes in 
peripheral relativistic heavy ion collisions. The coherent contribution, i.e, 
the excitation of a giant dipole resonance via meson exchange is shown to 
be negligible as well as interference between coulomb and nuclear excita- 
tion. Incoherent nucleon knock out contribution is also derived suggesting 
the nature of the nuclear interaction in this class of processes. We also 
justify the simple formulae used to fit the data of the E814 Collaboration. 



1 Introduction 



The important role played by the electromagnetic interaction in single nucleon removal 
processes induced by Relativistic Heavy Ions (RHI) has long been known. Dissociation 
of RHI by the Coulomb field of the target nucleus was first reported by Heckman and 
Lindstron [fij using for that purpose radiochemical techniques. Many other accelerator 
experiments have followed, using several combinations of heavy projectiles and targets 
with incoming energies E > 1. GeV/nucleon. The most conspicuous evidence of such 
electromagnetic processes is their large cross section and approximate dependence with 
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Z 2 , Z being the atomic number. Within this context it is very important to have a 
quantitative model which allows one to isolate and identify possible contributions arising 
from different processes due, for instance, to the strong interaction. 

Several simple models have been used to separate the nuclear contribution the most 
popular among experimental groups being the Limiting Fragmentation Model (LFM). 
This model assumes that the yield of the particular fragment from the target due to 
nuclear interactions will be independent of the beam, except for a geometrical factor 0- 
[[iJ . Based on this model a parametrization of the nuclear contribution is performed and 
the electromagnetic contribution extracted from the data. More recently in Brookhaven 
the same type of process is investigated || , Q . Their data analysis assumes the one nucleon 
removal to be basically due to Coulomb dissociation (parametrized as bZj*, b and c being 
free parameters) and to nuclear effects (parametrized as the sum of the heavy ions' radii 
(a (A^/ 3 + Ap 3 ), a being a free parameter). In this context an essential question arises: is 
it possible to give a sound theoretical basis for the contribution of the strong interaction in 
the one nucleon removal process? Can one theoretically justify the E814 prescription, for 
the nuclear cross section, quantitatively? Answering these questions are the main purpose 
of the present work. 

We consider two mechanisms which could be responsible for one nucleon removal me- 
diated by the nuclear field: the excitation of a giant dipole resonance via vector meson 
exchange and the direct incoherent knock-out contribution. In the first case, i.e., for the 



collective effects we use the framework of ref.[10] and derive analytical expressions for pure 
Coulomb (rather well known), pure nuclear and interference cross sections (here obtained 
for the first time). We show that the theoretical cross sections involving the strong inter- 
action are independent of the bombarding energy, depend very sensitively on the meson 
mass and display a different qualitative behavior as compared to the empirically predicted 
one from E814. Moreover for realistic values of the coupling constant and mesonic masses 
these cross sections are too small. We are therefore led to the conclusion that one nucleon 
removal through the excitation of a giant dipole resonance by means of the nuclear field 
cannot be the mechanism responsible for the observed nuclear contribution. 
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We next investigate the possibility that one nucleon is removed by means of and 
incoherent process (fragmentation reaction). The adequate theoretical tool for this purpose 
is Glauber's multiple scattering theory p| . It has long been known both for its simplicity 
and amazing predictive power. One can find copious examples in the literature where 
such theory allows for a simple physical interpretation of experimental results as well 
as their quantitative analysis ([|I^]-[|l4|]). We therefore work in such scheme. In fact, 
fragmentation reactions of the type discussed here have already been successfully analyzed 
in the framework of Glauber's theory: in one nuclear removal reactions, the momentum 
distribution of the outgoing fragment has been shown to reflect the momentum distribution 



of the nucleon which is removed from the surface of the projectile nucleus [14|. This is 
precisely the physics behind the incoherent nuclear one nucleon removal, as we shall show 
here. The target dependence so obtained corresponds precisely to the empirically fitted 
one. We are thus able to reproduce the results of the E814 experiment and to give a 
physical interpretation of the empirical fit. The calculation presented here, we believe, 
provides for a good theoretical prescription as to the quantitative nuclear contribution to 
the process. 

This paper is organized as follows: in section 2 the Coherent Effects, i.e, the excitation 
of a Giant Dipole Resonance via one photon and one vector meson exchange are analyzed. 
Analytical expressions for the pure nuclear and interference contributions are derived. 
Section 3, presents the study of the Nuclear Incoherent Effects, using Glauber's formalism. 
Finally we compare our results with E814 data showing very good agreement. 



2 Nuclear Coherent Effects 



Recently the E814 Collaboration || has reported on the total cross section for the 
excitations of the giant dipole resonance by means of peripheral heavy ion collisions. The 
reaction considered was 2S Si into p + 27 Al at an energy of 14.5 GeV/nucleon bombard- 
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ing different targets (C, Al,Cu, Sn, Pb). The experimental data suggest a non negligible 
amount of nuclear effects and show that the final-state energy, is peaked near the isovector 
giant dipole resonance in 28 Si. 

The pure Coulomb contribution is quantitatively the most important and has been 
previously obtained by several authors by different methods (Refs. |l(|, [p^j , []i"6f| ) . On the 
other hand the pure nuclear and interference cross sections are not as well understood (see 
ref. 0). 

In this section we study the excitation of a giant dipole resonance via one photon and 
one meson exchange and assume the one nucleon removal to be entirely due to this process. 
The study of vector and scalar meson exchanges in a covariant formalism has been 



performed in detail in ref. [18] for inelastic proton- nucleus scattering. In this reference it 
is shown that within the relativistic eikonal approximation one can sum an infinite series 
of Feynman diagrams and the resulting scattering amplitude is strongly dominated by the 
vertex where the excitation takes place. We therefore restrict ourselves in our derivations 
to the lowest order Feynman diagram contributions for the excitation of a giant dipole 
resonance via vector meson and photon exchange. 

We show that in the ultra-relativistic limit it is possible to derive very simple analytical 
formula for the three contributions (pure Coulomb, pure nuclear and interference), to 
understand the physical reason why the pure nuclear contribution is so much smaller than 
the Coulomb one and to estimate the contribution of interference. 

2.1 Collective Coulomb plus Nuclear Excitation Cross Section 



Formally, the first order contribution to collective nuclear excitation via one photon 
and one meson exchange between projectile and target nucleus is given in the Center-of- 
Mass System (CM.) by JTH]: 
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dfl 
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(1) 



where a is the fine-structure-constant, is the total center of mass energy, the Z's and 
M's are atomic and mass numbers of each nucleus, respectively, q is the four-momentum 
transfer and q A (<7b) its spatial part in the proper system of the nucleus A (B) (ref. |To|l ). 
g v is the strong coupling constant and fi is the vector meson's mass. 

The F l A are the Fourier transforms of transition current matrix elements, given by: 



F l A (q A ;X A ,m A ) 



C 



d A r < E A ;X A ,m A \j l A (r)\Ag.s > e 



< E* A ; X A ,m A \j A (q A )\A 



9-s 



> 



(2) 



with C = eZ A for the Coulomb contribution and C = g v for the strong contribution. We 
can write the similar expression for the nucleus B. 

In the above equations, currents and nuclear states are represented in their proper sys- 
tems, and \E; A, m > stands for the final state of nucleus A or B with angular momentum 
A whose projection in the incident beam direction (e' z ) is m. J2f means the summation 
over final states. The Lorentz transformation matrix is given by 
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(3) 



where (3 is the velocity of nucleus A in the rest frame B, and 7 is the associated Lorentz 



factor (7 = y/l -P 2 ). 

In order to evaluate the current matrix elements, it is convenient to use the coordinate 
system whose z-direction is taken to be the direction of the momentum transfer seen by 
the nucleus. The cartesian base vectors {e x ,e y ,e z }, where z = q/\q\, are related to the 
original ones {e' x ,e' y ,e' z } by 
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(4) 



where 6 and 93 are polar angles of z with respect to the base vectors {x',y',z'}. The 
corresponding spherical base vectors {eli, eb, e+i} and {e' +1 , e' , e'_ ± } are related to each 
other by 



e+i ^ 






eo 




e'o 


e_i y 







(5) 



where (99, 6*, —99) is the usual rotation matrix for total angular momentum j = 1. Using 
these relations we can express the four-vector current in terms of spherical components as 



(6) 



where the 3x3 matrix U transforms spherical vectors into cartesian ones. The convenience 
of the spherical basis now becomes apparent: the selection rules for the current matrix 
elements are 
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< E*;\,v\p\0,0 >=< p > S vfl , 

< E*;\,v\j +1 \0,0>=<j +1 > S Vt+ i, 

< E*;\,v\j \0,0 >=< j > 6 ufi , 

< E*; A, v\j-!\0, >=< j-! > S v -! =< j+i > S v -u 
where v refers to the projection of angular momentum in the three q directions. 



(7) 



6 



We can now write the matrix element in eq. (1) as 
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q (the argument of D(q)), symbolically represents the Euler angles associated to the rota- 
tion required to transform e' z to e z in each nucleus. 

In the case where the nucleus B stays unexcited, we obtain 



da 
dtt 
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= 2 Fc ' <PB> I V 
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m=0 



F c \ <p B > | V 



cos 6 A {< PA > -/5cos 9 A < jo >} - -p sin 2 9 A < jf > | 2 (10) 



where 9a represents the polar angle of q\ m the rest frame of A and similar expressions 
can be constructed for the strong and interference part with the respective substitution of 
F c by F s and F cs . 

The next problem one has to face in the present formulation is the inclusion of nuclear 
absorption effects. The way to circumvent this problem is the following: we first write the 
total cross-section in terms of an integral over the transverse momentum transfer, 
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where -JF~ C = 2qZ^Z^ and y/F~ s = g v . We then introduce a Fourier transform 



where M(q±) denotes 
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and p stands for the projectile's center-of-mass momentum. The total cross section can 
now be written as 



J d%Tr[f\b)f{b)} 



(14) 



The above formula can be interpreted as the impact parameter representation of the cross 
section. It is now possible to estimate nuclear absorption effects in a simple way: we 
attribute those parts of the integral coming from all values of b smaller than the sum of 
radii of projectile and target to the absorption coming from nuclear interactions which do 
not contribute to excite the collective mode. Thus, with b m i n = ro(Ap 3 + Ap 3 ), we get: 
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lta 2 (<u)-f 2 )(<! 2 W)-;< 2 ). 



Tr 



q±)e 



(17) 



and 



a. 



poo r r 

/ d 2 b / d 2 q ± / d 2 q' ± 



Tr 



(18) 



correspond to the Coulomb excitation process, to the pure strong contribution and to 
the interference between two, respectively. 



2.2 Results 



In order to obtain analytical results we now need a model for the matrix elements 
corresponding to the excitation of a giant dipole resonance. We use a macroscopic model. 
T. Suzuki and D. J. Rowe (T^j have derive a sum rule for current matrix elements which 
corresponds to the incompressible fluid model. In it, the giant resonance state is repre- 
sented by a unique level of given multipolarity and the transition matrix elements of the 
current from the ground state to this level is assumed to be of the form 

< g.s.\J{r)\E*,\v >= ^Po(r) V {r X Y Xu } (19) 

where po(r) is the ground state density distribution, m the nucleon mass and N the 
normalization constant related to the sum rule value S as, 

N = ^E*/S (20) 
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with 



S = SXAR 2X ~ 2 /87rm 



(21) 



The Fourier transform of eq. (19) gives the required matrix element. After some algebra 
we obtain 



<p>~- 



AV4vr(2A + 1)A / drr^poi^jx^iqr) 



2mi x E*/\q\ 

where j\+i(qr) are spherical Bessel functions and 
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(24) 



In this model, magnetic current is neglected. It is known that for giant resonances, the 
contribution of magnetic current is small. For very high energies (7 >> 1) and in the low 
momentum transfer limit {qR « 1), we may approximate the nuclear matrix element in 
eq. (22) as 



<p>~- 



AE* \q A \ 



2m E* 

Furthermore, the four-momentum transfer squared is given by 



(25) 



In this limit, 



q 2 ~ 



* \ 2 



E*\ 
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and 
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MAqA ~ -(AE*/2m) 1/2 q ± e- 
P 



(27) 
(28) 
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and it is possible to carry out the integrals in q±_ and q' ± as well as over the impact 
parameter b. We obtain for the Coulomb contribution 



a c = 3.5 x 1Q~ 3 Z% Na ^ In [6/£ + 1] mb. (29) 

where £ = E*b min /(3~f and 5 = 0.681085 .... 

The pure strong interaction and interference contributions can be obtained in an anal- 
ogous manner. We get 

cr s = 26.65(27r 2 )^^fV 2 ^""™ mb (30) 

and 



Pes = — Z B A^ mb (31) 

10m \/^b Tl 



As discussed before, expression (29) has been obtained before by several authors by dif- 
ferent methods. The results for the strong (a s ) and interferency part (a cs ) , eqs. (30) and 
(31), we believe to be novel. 

From the qualitative point of view it is clear that the theoretical prediction for the 

i i 

bmin = r o{A\ + A^) dependence of the total strong cross section (eq.(30)) is in contra- 
diction with the empirical parametrizations obtained in the E814 collaboration given by 
(Ref.§): 

a exp . s = a{A l f + A)( 3 ) (32) 
where a = 1.34±0.19/m 2 . The theoretical cross sections depend very sensitively on the re- 

2 

alistic values of the exchanged meson mass ^i(MeV) and of the coupling constant = 0.9 
I ffifl . For these reasons the nuclear and interference contributions to the collective excita- 
tion become negligible. Also as expected these two contributions are energy independent. 
The conclusion then is that the collective effects are essentially due to coulomb excitation 
and that the proposed mechanism is in rather good agreement with the data (see table 
1, columns 2 and 4). Notice, however, that the experimental values for the strong part 
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target 


0~exp— c 


0~exp— s 


0"c 


&CS 


OS 


C 


6.81 ± 1.03 


7.13 ± 1.00 


6.34 


1.13 x 10" 


-2 


7.00 x 10" 


-3 


Al 


28.15 ±4.88 


8.08 ± 1.14 


30.71 


6.68 x 10" 


-3 


1.70 x 10" 


-3 


Cu 


122.79 ±24.16 


9.40 ± 1.33 


131.04 


2.63 x 10' 


-3 


1.69 x 10- 


-4 


Sn 


333.87 ± 70.92 


10.64 ± 1.50 


366.79 


9.14 x 10" 


-4 


1.51 x 10- 


-5 


Pb 


828.25 ± 187.65 


12.00 ± 1.70 


927.11 


2.50 x 10- 


-4 


9.01 x 10- 


-7 



Table 1: Experimental Coulomb (cr exp _ c ), Experimental Strong (<J exv - s ), our Coulomb 
(a c , eq.29) interference (a cs , eq.31), and our Strong cross section in the last column (a s , 
eq.30), with g v = 3.36 and \i = 300Mey. Units in mb. 

of the interactions are non negligible. This indicates that one should rather investigate 
incoherent contributions to the nuclear part. In fact it can be shown that the incoherent 
knock-out process gives both qualitative and quantitative agreement with the empirical 
findings. 

2.3 Conclusions for Coherent Effects 

In this section we use a covariant formulation to describe the excitation of a collective 
mode via photon and meson exchange. Besides the well known total Coulomb contribu- 
tion, analytical results are obtained for the total strong interaction contribution and the 
interference cross section. We show that the fact that the strong interaction is mediated 
by massive mesons drastically reduces its contribution to the excitation of the giant reso- 
nance. Both the pure strong and interference contributions are independent of the incident 
energy. 

Moreover the theoretical prediction for the strong part of the collective total cross 
section is in qualitative disagreement with the experimental findings. 

Therefore we conclude, that the electromagnetic effect is the only responsible for the 
processes where coherent effects are produced, as a giant dipole resonance. The strong 
interaction contribution is coming from incoherent effects such as one nucleon knock-out. 
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This point is investigated in the next section. 



3 Nuclear Incoherent Effects 



In this section we analyse Glauber's formulation of the one nucleon removal process 



according to reference [14], and give a theoretical expression for the total nuclear cross 
section. We show its compatibility with the empirical formula used in the analysis of the 
E814 experiment ||. 

The physical process considered is a nucleon knock-out off the projectile within the 
context of Glauber's formalism. 

3.1 The Nuclear One Nucleon Removal Cross Section 



In this section the contribution of the strong interaction for the one nucleon removal 
process is calculated using Glauber's multiple scattering formalism [14|. The notation is 



summarized in the following equation, valid in the projectile's rest frame: 

A Z{\0,^ >} + T{\ -p ,9 >} -V 4 - 1 Z{\k,4> a >} 

+n{\p, np >} + T'{\ -p - q, 9fs\ >} (33) 

Before the reaction, the projectile A Z, with intrinsic wave function is at rest and the 
target approaches with momentum — p~Q. After the interaction, a fragment with A — 1 
nucleons is detected in a particle stable state (ft a with momentum k. The nucleon with 
momentum p and scattering wave function rjp, and the final state of the target, are 



unobserved. The cross-section corresponding to the equation (33) is given by [13| 
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I fp e -i^/ :@p]1 _ ft (i-r y )|v*;e )| a (34) 

In the above formula Yij{xi + b — yj) are the profile functions for collisions between a 
projectile nucleon i located at X{ and a target nucleon j at yj; ifj and yj refer to the 
respective centers of mass and the impact parameter b is the transverse distance between 
projectile and target. 

The fragmentation cross-section corresponding to the incoherent process is derived in 
detail in reference We quote the result 

^ = J D(s)d 2 s J dz J d 3 kW{s,z,k) (35) 

where s is the perpendicular distance of each nucleon to the center of the nuclei. The 
distorting function D(s) contains the reaction dynamics 

£)(g) — J ^2^ e -27mxFT(fe)-2/mxT(?+s)^ e o-iVJv(0)T T (fe+s) _ -g 

iXFT(b) ~ -{<t>a®a\ r mn|<^a@o) (37) 
m>l 



where 



(b + S) « -(G | £ r ln (6 + s - 4)|G ) (38) 



T T (B) = A T J dzp(B,z) (39) 
with p(B,z) is the nucleon density function, and 



is the nucleon-nucleon cross section. 



d 2 q ^r b djN (4Q) 

erg 
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As discussed in reference (||14||) shown below, D(s) localizes the reaction to the nuclear 
surface. The Wigner transform W(s, z, k) contains the momentum distribution of the 



removed nucleon. 



In the calculations we have used Gaussian densities 



(41) 



where for light nuclei we have adjusted ao in order to reproduce the root mean square 
radius. For heavier nuclei (A > 40) N and ao have been adjusted in order to reproduce 
the tail of the Woods-Saxon density. 

The theoretical results to be shown are based on the result 



kmin corresponds to the Coulomb barrier cut-off parameter, in the case where the nucleon 
removed is charged (<7tv is the nuclear cross section). 

3.2 Results - E814 Brookhaven Experiment 

In the E814 Experiment was reported the results of the EMD of 28 Si into p+ 2 Al, and 
Pb, Sn, Cu, Al and C were used as targets. The total cross section for producing p + 27 Al 
was measured. Presenting the first measurements of the final-state energy spectrum (for 
all targets) which is peaked near the isovector giant-dipole resonance in 28 Si. The incident 
energy was 14.5 GeV/nucleon. A detailed description of the experiment can be found in 
ref.||. The empirical fitting was for the electromagnetic part: 




(42) 



with 




(43) 



exp—c — 



bZ^ 



(44) 
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and for the nuclear part: 

a exp _ s = a(4 /3 + A 1 / 3 ) (45) 

with a = 1.34 ± 0.19/m 2 ., b = 0.23 ± 0.05/m 2 . and c = 1.8 ± 0.06.. Using this fitting the 
authors showing a good agreement for all ions. We calculate the nuclear incoherent cross 
section (eq.42), using the corresponding cut-off for the Coulomb barrier. 

Before a comparison between the theoretical results and the E814 data a word about 
the adequacy of such theoretical description to this specific set of data is in order: in the 
E814 experiment events which correspond to an energy transfer above 200 MeV have been 
excluded precisely in order to avoid contributions coming from the nuclear interaction. 
In Glauber's model the energy transfer is realized by means of nucleon-nucleon collisions 
which are assumed to occur essentially in the impact parameter plane. The average mo- 
mentum transfer is 400 MeV jc which corresponds to and energy of about 80 MeV . Also 
the highest probabilities occur for smaller momentum transfers. We can therefore not 
exclude such processes from the experiment. The qualitative and quantitative agreement 
with experiment is, in the end what strengthened our conclusion. 

In Table 2 we can see the nuclear cross section fitted by the experimentalists and our 
theoretical calculation of the nuclear part. We find good agreement. In this way our 
conclusion is that the processes is due to incoherent effects. Our numerical result based in 
eq. (42) gives a nuclear cross section proportional to the sum of the radii of the target and 
projectile, this is showing in the third coulumn of this table where the ratio ^1/3^1/3 ^ ^ s 
shown to remain constant. 

We demonstrated that in this class of process, the interference cross section between 
Coulomb and strong interaction is negligible. The nuclear part is due to an incoherent 
process. In the context of Glauber's theory, we show that the nuclear cross section is 
proportional to the sum of the radii of nuclei, in perfect agreement with the experimental 
fitting used in E814 Collaboration. 

In table 3 we can see the total cross section obtained in the E814 Collaboration in 
comparison with our theoretical results for the electromagnetic coherent and nuclear in- 
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A T 


a exp - s (E8U) 


a N 


a^/(4 /3 + 4 /3 ) 




7.11 ± 1.00 


8.11 


1.5 


27 Al 


8.04 ± 1.14 


9.24 


1.5 




9.34 ± 1.33 


11.12 


1.5 


U8 Sn 


10.59 ± 1.50 


13.04 


1.5 


208p b 


11.95 ± 1.70 


15.23 


1.6 



Table 2: First column shows the parametrization used in E814 (nuclear part), second 
column our results, third column shows the proportionality with the sum of the radii. 



A T 


a exp _ t (E8U) 


0~N + °- cs 




13.92 ± 2.03 


14.45 


27 Al 


36.19 ±6.02 


39.95 


63 Cu 


132.13 ±25.49 


142.16 


118 Sn 


344.46 ± 72.47 


349.83 


208p b 


840.20 ± 189.35 


942.34 



Table 3: First column shows the total cross section data of the E8I4 Collaboration, second 
column show the sum of our results for the nuclear and electromagnetic cross section. 



coherent parts. Good agreement is found. 



4 Conclusions 



In the present work the nuclear and interference nuclear-coulomb contribution to one 
nucleon removal processes in peripheral RHI collisions are studied. Firstly its coherent 
contribution is analyzed assuming that a giant dipole resonance is excited via one vector 
meson exchange. To lowest order simple analytical expression for the electromagnetic, 
nuclear and coulomb-nuclear interference contributions are derived. 

For the electromagnetic part a cross section proportional to « Z 2 is obtained, as is 
expected experimentally and theoretically. The interference and coherent nuclear part 
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are small and in qualitative disagreement with the data. We therefore conclude that the 
nuclear contribution in experiment of one nucleon dissociation, results from incoherent 
nucleon-nucleon collisions. Therefore, in section 3, the incoherent one nucleon knock-out 
is studied within the context of Glauber's formalism. The total nuclear contribution so 
obtained is shown to be proportional to the sum of the radii of the nuclei as could be 
intuitively expected. Quantitatively speaking the incoherent process is shown to be of the 
order expected in experiments. 
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